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Abstract 

Measurements of various physical quantities in the symmetry broken phase of 
the one component lattice (/>| with standard action, are shown to be consistent with 
the critical behavior obtained by renormalization group analyses. This is in contrast 
to recent conclusions by another group, who further claim that the unconventional 
scaling behavior they observe, when extended to the complete Higgs sector of the 
Standard Model, would alter the conventional triviality bound on the mass of the 
Higgs. 



*On leave from Eotvos University, HAS Research Group, Budapest, Hungary 



1 Introduction 



The continuum limit of the lattice regularized (j)'^ theory in 4 dimensions with n > 1 
components is thought to be trivial. The evidence comes from renormalization group 
studies and numerical simulations; remarkably however there is still no rigorous 
proof. 

After a period of intense numerical investigations about 15 years ago, interest 
in this model within the lattice community dwindled. Occasionally however some 
authors have questioned the usual procedures of renormalization in the broken phase. 
In their latest paper [J Cea, Consoli and Cosmai claim to present further numerical 
evidence for their scenario. From the presentation of the results in their paper a 
neutral reader would, accepting their analysis, indeed conclude that the conventional 
wisdom (CW) concerning the critical behavior in the lattice (j)^ was incorrect. The 
purpose of this paper is to point out the deficiencies of their analyses, and to show 
that the numerical data is completely consistent with the conventional picture. 

2 The lattice theory and renormahzation group predictions 

The standard lattice (j)^ model is characterized by two bare parameters k, A; the 
action is ^ 



There are two phases separated by a line of critical points k = Kc(A). For k > Kc(A) 
the symmetry is spontaneously broken and the bare field (p{x) has a non-vanishing 
expectation value v. As usual this is defined by first taking the thermodynamic limit 
in the presence of an external uniform magnetic field h and then letting this field 
tend to zero 

V = lim lim (0(0))/, . (2.2) 

This study will restrict attention to the vacuum expectation value and the 
connected two-point function 



A renormalized mass and a field wave function renormalization constant are 
defined through the behavior of Fourier transform of G{x) for small momenta: 
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s = J2 -^^Yl Hx)<P{x + A) + Hxf + K^ixf - if 



(2.1) 



X f_l=l 




(2.3) 



G{k) 



(2.4) 



^We use the notations in refs. i2ii3i. 
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The susceptibility x ^-iid the second moment ^ are defined through 

X = Y.'^{x) = Zjml, (2.5) 

X 

H = x^G{x) = SZj,/mi . (2.6) 

X 

We also define the normalization constant associated with the canonical bare field 
through ^ 

= 2kZ^ = 2Kmlx . (2.7) 

In the framework of perturbation theory correlation functions of the multiplicatively 
renormalized field 

4>^ {x) = Z-^'^<p{x) , (2.8) 

have at all orders finite continuum limits after mass and coupling renormalization 
are taken into account. Correspondingly a renormalized vacuum expectation value 
is defined through 



v^ = vZ^^/\ (2.9) 

Here we adopt the generally accepted assumption that the structurally same mul- 
tiplicative renormalization is required to define the continuum limit of the theory 
non-perturbatively. We see no evidence for the claim in ref. ^ that the vacuum 
expectation value and the fluctuating part of the bare field should be renormalized 
differently. Finally a particular renormalized coupling is defined by 

g^=3mi/vl=3mix/v\ (2-10) 

as is one popular choice in the perturbative framework. 

The commonly accepted picture of the critical behavior of the theory is mainly 
due to the work of the Saclay group [1]. In their framework the renormalization 
group equations predict that the mass and vacuum expectation value go to zero 
according to 

ma OC t1/2[1j^(^)[-1/6^ (2.11) 
Vn OC Ti/2[ln(r)li/3, (2.12) 

for T = k/kc — 1 — > 0, and correspondingly the renormalized coupling is predicted to 
go to zero logarithmically which is the expression of triviality. The critical behavior 

^In various papers a different notation is employed and the quantity 16kx^/m is denoted by Zn. 



2 



in the broken phase is conveniently expressed in terms of three integration constants 
C- ,i = 1, 2, 3 appearing in the critical behaviors: 

m^. = Ci(/3i5Ry'/'Vi/^^5Ml + 0(<7R)} , /5i = y|^> (2-13) 
= C'.{^-l,^ + Oi,l)}, (2.14) 

K-K, = ^C>|5r'^'{1 + 0(5r)} . (2.15) 

In ref. jHj these constants were estimated by relating them to the corresponding 
constants Cj in the symmetric phase. These in turn were computed by integrating 
the renormalization group equations with initial data on the line k = 0.95kc{\) 
obtained from high temperature expansions. 



3 Some proposed criticisms of CW 

In ref. ^ the authors present two quantitative arguments against the conventional 
wisdom. Their work considers just the Ising limit (A = oo). Firstly they claim 
grows logarithmically as one approaches the critical point instead of going to a con- 
stant as predicted by the RG (see Eq. H2.14() ). Unfortunately they only measure v, x 
which is not sufficient to determine Z^. To overcome this deficiency they compute 
Zji = Jn-fnputXi where x is their MC measurement at a given value of k but minput 
is an estimate of at the same k taken from Table 3 of ref. j3I (referred to as T3 
below). This "composed" quantity indeed seems to grow as the critical point 
is approached. Accepting that their measurements of x correct (and indeed we 
agree with their values on the lattices we checked), the crucial question is whether 
the estimates m-input of are reliable. 

At this stage two details of the analysis in ref. [3] must be appreciated. Firstly 
the values of k cited in the last column of T3 are written using Kc obtained from the 
same procedure of analysis for all values of A ^. For the particular case of the Ising 
model there are (presumingly) more accurate values; some comparisons are made in 
Tabled Hence, even if one accepts the estimated errors in T3, estimates of for 
a given k obtained using this table naively are subject to further large uncertainties. 

Secondly it is probable that the systematic errors in T3 have been underes- 
timated. This stems from a probable underestimate on the size of the systematic 
errors on the cited values of the integration constants C'^; e.g. the estimated values 
of the constants obtained by using 2-loop or 3-loop expressions in the RG equa- 
tions differ considerably. This question was briefly addressed in the 3rd paragraph 

^In hindsight it would have been better to have tabulated estimated values of k — Kc 
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ref. 



0.07475(7) 

0.074834(15) 

0.074848(2) 

0.074851(8) 

0.07487(5) 



(21 

m 
m 
m 



Table 1: Estimates of Kc in the 4d Ising model. 



sect. 5.2 of 3 and later stressed by Peter Hasenfratz [S] ^. 

The simple outcome of the discussion above is that to compute and 
independently one needs accurate measurements of both x, A*- In the next sections we 
describe such measurements and find that the computed values of are consistent 
with the RG expectation in the (j)^ theory. For the Ising limit the measured values 
of rriR are considerably lower than the corresponding estimates minput in 

The second presented "evidence" against CW in ^1_ concerns the quantity v'^x- 
Defining 



in the RG framework v'^x behaves as oc £ as one approaches Kc, whereas in their 
scenario they expect a more singular behavior oc i'^. Figure 2 in gives the im- 
pression that the data strongly supports their prediction. However, this impression 
is completely false! The figure is obtained by 2-parameter fits of v'^x with functions 
Aln^{n — Kc) with p = 1,2 and Hq a fit parameter. But such fits are of course mean- 
ingless without including a scale of the logarithm. Once this is done (as described 
in subsect. 4.4), one gets a beautiful fit to the data also for the function expected 
from the RG group. 

4 Ising MC simulation and results 

In this section we present our simulations of the 4d Ising model on a hypercubic lat- 
tice with periodic boundary conditions in all directions. We consider only geometries 
with volumes V = L^T. 

^This critique was discussed in more detail for the physically relevant case of 0(4) where 
fortunately the various loop estimates are closer. 



ln{n — He) 



(3.1) 
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4.1 Simulation algorithm 



For updating the configurations we used the Swendsen-Wang cluster method [TT] . 
Beyond practically eliminating the problem of critical slowing down, the method 
allows to use improved estimators which reduce significantly the errors of the mea- 
sured quantities. Since the application of the cluster improved estimators have some 
specific features for the broken phase, we discuss briefly here the applied procedure. 

The first step in the SW cluster algorithm is to put bonds between neighboring 
spins of equal signs with probability p = 1 — exp(— 4k), while no bonds are put 
between spins with opposite signs. In the second step one identifies the clusters, 
the set of spins connected by bonds. Obviously, spins within a cluster have the 
same sign. Denote the number of clusters by v, and the number of sites in the 
z-th cluster by n^, i = 1,... ,z/. Assume for convenience that the numbering is 
chosen so that rii > n2 > ... > n,^. In the updating step one flips all spins in 
the i-th cluster with probability 1/2. All resulting 2'^ configurations appear in the 
equilibrium distribution with equal probability. The cluster improved estimator 
replaces a quantity by its average over these 2'^ configurations. In particular the 
product 4){x)4){y) is replaced by +1 if x and y belong to the same cluster, and by 
if the two sites belong to different clusters. 

4.2 Determination of v and x 

The difference between the symmetric and broken phase shows up in the distribution 
of the cluster size. In the symmetric phase (at finite correlation length ^) the typical 
size of clusters is determined by ^ and remains finite and independent of L for 
S> ^. This is obvious from considering the spin-spin correlator G{x) at large 
distances |x| since this number is equal to the probability that the two sites are 
within the same cluster. As a consequence, this probability goes down exponentially 
with the distance. It is easy to show that in the symmetric phase the susceptibility 



This involves that the size of the largest cluster grows proportionally to the volume, 
{rii) oc V while the distribution of n^, i = 2, 3, . . . is not affected by the volume for 
L,T ^ (and is characterized by the correlation length). In this case (for V — > oo) 
the vacuum expectation value is given by 



is given by x = Ei ^i) ■ 



In the broken phase (in infinite volume) {(f){x)4>{y)) — v 



for \x — y 



oo. 



V = 




(4.1) 



and the susceptibility by 




(4.2) 
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In sufficiently large volumes the definition Eq. H4.1|) of the vev coincides with the 
conventional definition Eq. ()2.2|) . For sufficiently large V in the broken phase one 
can choose an external field h such that h{ni) = hvV ^ 1 and at the same time 
hn2 <^ 1. Therefore in the presence of the magnetic field the spins in the largest 
cluster are frozen to value +1 while h has practically no influence on the orientation 
of the other clusters. Eq. ()4.1() provides a convenient definition of in a finite 
volume. For large enough volumes this definition also coincides practically with the 
Binder's definition jl3j using the absolute value of the total magnetization: 



"^Binder 




which is the definition employed in ref. ^ . For a thorough discussion of various finite 
volume definitions of the vev (and more general finite volume effects in this model) 
we refer the reader to ref. Assigning random signs to the clusters i = 2, . . . ,v 
the sum of spins in these clusters has a variance X]r=2 '^f ~ '^X (cf- Eq. (|4.2j) ). 
while the value of the spin in the largest cluster is ni ~ vV . The probability that 
the random component has the same magnitude as the constant component ni is 
suppressed at least by the factor exp(— Fv/x^) . 

Note that for the case oiV = L^T, L = fixed, T oo the situation is somewhat 
different. In this quasi-one-dimensional case the large clusters have finite length 
in the i-direction, although much larger than L. Their typical length is given by 
the inverse of the energy gap between the lowest even and odd eigenstates of the 
transfer matrix Here we shall not discuss this geometry, note only the obvious 
fact that any quantity defined through spin correlators can be expressed in terms of 
expectation values of the cluster sizes rij. 

Apart from v and x we also measured the time-slice correlation functions 

5(t) = j;G((t,x)), (4.4) 

X 

to compute the exponential mass and the Fourier transform G{k) at momenta k 
along the t-axis. 

Typically we performed more than 400K sweeps at each k point. The autocor- 
relation times for the vev and the two point function S{t = T/4) were monitored 
and found reasonably small e.g. for the lattice with the largest correlation length 
which we measured («; = 0.0751) they were ~ 8, ~ 3 respectively. 

4.3 Hacw data and analysis 

To illustrate the quality of the raw data, and for later reference, in Table I^J the 
time-slice correlation function S{t) is given for V = 48^ at k = 0.0751. 
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17.788(31) 


7 


5.320(27) 


14 


1.666(25) 


21 


0.646(29) 
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14.913(31) 


8 


4.492(26) 


15 


1.421(25) 


22 


0.600(29) 
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12.529(30) 


9 


3.796(26) 


16 


1.217(26) 


23 


0.573(30) 


3 


10.540(30) 


10 


3.211(25) 


17 


1.048(26) 


24 


0.563(30) 


4 


8.875(29) 


11 


2.720(25) 


18 


0.909(27) 






5 


7.478(28) 


12 


2.306(25) 


19 


0.798(27) 






6 


6.306(28) 


13 


1.958(25) 


20 


0.711(28) 







Table 2: The (subtracted) time-slice correlation function S{t) for the Ising case at k = 0.0751 
on a 48'' lattice. Note that the quoted errors for this quantity are much smaller than that for the 
subtracted term «^L^ = 2711.1(1.2). 

In Table 01 we collect data of the values of v,x for various values of k,L from 
two papers jl2lll together with our own. The simulation at the largest value of 
K = kq was performed to test our programs and analyses by comparison at a point 
where the low temperature expansion is considered quantitatively reliable [S]. The 
agreement between the data and analytic computations is indeed satisfactory. As for 
the values of k closer to Kc, firstly one observes excellent agreement of the raw data 
on those lattices measured by different groups. Thus we are confident that the raw 
data (albeit sometimes obtained by different methods) can be trusted. Secondly at 
values of /t where different volumes were measured there are no signs of significant 
finite volume effects. 

The second moment mass can be determined directly from the Fourier transform 
of S{t) (i.e. G{k) with momenta along the t-axis) using Eq. (|2.4() which involves a 
determination of the slope from the available discrete values of A;. To avoid possible 
discretization errors due to the finiteness of the lattice size T we also used the 
following procedure to determine and Z^. We first computed the exponential 
mass m characterizing the exponential fall-off of the time-slice correlation functions 
at large separations, using one-mass and constrained two- mass fits (where the second 
mass was required to satisfy ?7i2 > 2m). For the one- mass fit we included only 
distances t with mt > 1. The results of the two fits were completely consistent, the 
central value of the 2-mass fit being slightly lower but the estimated errors larger 
than those of the 1-mass fit. In Table lU we quote only the outcome of the 1-mass 
fit. Note that the value of rriL is in each case large > 8, and hence finite volume 
effects are expected to be smaller than our statistical errors. An estimate of the 
infinite volume second moment mass was then computed from the sums Ylt "^(O 
and ^tt^S{t), where for |i| < Tq (some large Tq < T/2) we took the measured 
correlations and for Tq < |t| < oo we estimated S{t) by assuming it has the form 
of a free lattice correlator with the previously determined mass m and amplitude. 
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The resulting estimates of from the two fits hardly differed; in Table 0] we quote 
the results from the 1-mass fit. The measurements of from G{k)~^ were also 
consistent with the values in the table (although statistical errors were typically 
a factor 2 larger), again indicating small finite volume effects. We also remark 
that the relations between the quoted and m are consistent with estimates from 
renormalized perturbation theory j3ll2j . 

At this stage we reach the goal of determining through Eq. (|2.5j) . The 
corresponding values of are tabulated in Table [l] For all values of k we find 
Z^ < 1 and consistent with the CW of tending to a non-zero constant as k ^ Kc- 
There is no signal of a logarithmic increase as favored by the authors of ref. pj 

What does not agree so well in our measurements with T3 is the value of g^i 
for a given m^. As mentioned before, the estimate of in ^ is sensitive to the 
integration constant lnC( = InCi + 1/6; the estimates in T3 are obtained with 3- 
loop RG equations and an input value of InCi = 1.5(2). The present measurements 
of 5r are better reproduced with a central value of ~ 1.2 for InCi. 

4.4 The quantity v'^x 

In Fig. 1 we plot fits of the data for v'^x with functions ^ 

Ailnf{K-K^)+A2, p=l,2. (4.5) 

To make our point, for the fits we just used the data of ref. Our data (the 
triangles) fall nicely on the fits and illustrate again the consistency of the data. 
Both fits are of good quality with chi^/dof ~ 1. The immediate conclusion is that 
the data cannot distinguish between the critical behaviors; to accomplish this one 
would have to get much closer to the critical point and treat much larger correlation 
lengths. 

Apart from the rather solid theoretical foundation, an additional strength of the 
CW is that the coefficient of the log (^i) is related to C2, which has been estimated 
from data in the symmetric phase. In fact conventional wisdom predicts a form: 



v'^X = ai 



25^ 
— m^ 
27 



+ «2 + O ( i ) , (4.6) 



with 



02 = ai [In + 2 In C[ + K] , (4.8) 
K = iln(^)-l(7 + 21n2) = -1.6317... (4.9) 



^just adding a constant in the case p = 2 is of course completely ad hoc, but the authors of this 
paper are actually only interested in the case p = 1 
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Using the values of the in one gets ai = 1.20(3) ,02 = —1.6(6). Fitting the 
entire data set with k < 0.0759 to the function (of the RG expected form) 



A, 



27 



+ A2 



(4.10) 



we obtain a good fit with chi^dof = 0.9 and Ai = 1.267(14) , A2 = -2.89(8) , Kc = 
0.074833(17). The fitted value of Kc is consistent with the best values given in 
Table [H Further the Ai are in reasonable agreement with the predicted values 
above. 



5 Field theory simulations at constant physics 

We have also studied the cutoff dependence of the renormalization constant 
along a line of constant IR physics, keeping the renormalized coupling g^^ (defined 
in Eq. ()2.10|) ) fixed as the renormalized mass in lattice units is varied. This is 
of course the behavior of interest in deciding how limits on the appearance of new 
UV physics follows from the triviality of the lattice model in the continuum limit. 
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This has been done at three points A = 0.3,0.6 and 1.0 (the Ising limit) ^ along the 
RG curve corresponding to (^r ~ 20 (cf Fig. 2 of The simulations in the field 
theory cases (A = 0.3, 0.6) were done with conventional Metropolis update code on 
lattices of size 64^, while the Ising point results were obtained using a cluster code 
on lattices of size 48^. 

The renormalization constant and the zero momentum renormalized mass 
were extracted from a fit of the inverse lattice propagator in momentum space 
for small lattice momenta: 

Gik)-^ = Z-\ml+P + 0{k^)). (5.1) 

In the cases (A = 0.3, 0.6), we have computed G{k) = (|(/)(/c)p) from the FFT (fast- 
Fourier-transform) (j){k) of the coordinate space field (/'(x). The FFT is performed 
every 40 Monte-Carlo sweeps: we find that the nonzero momentum modes of G{k), 
for which the vev is irrelevant, have an autocorrelation time ranging from 50 to a few 
hundred sweeps. For A=0.3 we have collected propagators for a total of 50K sweeps 
(1250 propagators), while for A =0.6 we have lOOK sweeps (2500 propagators). 

Only the lowest 4^-1=255 modes are included in the fit: the mode corresponding 
to zero momentum is omitted. We then perform an uncorrelated (diagonal chi- 
square) fit of the 34 data points (corresponding to different values of fc^) to Eq. (|5.1() . 
The chi-squared/degree of freedom for the A =0.3 (resp. 0.6) point was 40/32 (resp. 
54/32). The results for Z^,m^ obtained by this procedure are indicated in TablelS 
Also shown is the field vacuum expectation value (unrenormalized) {(j)), from which 
a renormalized coupling y^, also shown in the Table, can be obtained. Finally, we 
give the susceptibility x ^^'^ the value 2Km|x which is the value of calculated 
from the zero momentum propagator, Eq. (|2.5jl . 

It is apparent from the results in Table El that the cutoff dependence of Z^ is 
extremely mild over the entire range covering a factor of 3 in cutoff holding the 
physical mass fixed, perfectly in accord with conventional renormalization wisdom. 
Note that the values of 2Km\x agree with the latter determination of within the 
errors, as it should be since the Fourier transform of the (subtracted) correlator, 
G{k) is a continuous function. This again shows that the determination of the bare 
vev is correct. 

As there may be some question of the extent to which the inverse momentum 
space propagator is adequately fit by Eq. 1)5. 1|) with two parameters, we have dis- 
played the quality of the fits for A=0.3, 0.6 in Figures 2, 3. There is certainly no 
evidence for any curvature up to the maximum momenta included in the fit, and 
indeed the chi-square/degree of freedom for the two-parameter fit is perfectly fine. 
In Table [51 for the Ising case we quote numbers for mR and Z^, obtained by fitting 

®for the definition of A see 
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Figure 2: [2KG(fc)]"^ at A=0.3, ~ 20. The triangles are the data from 50K sweeps 1250 
propagators. The fit is with Zi{,,mR given in Table El 



the first three k ^ values (with k along the time-axis) of the inverse propagator. 
These numbers differ only slightly from those given (for the same k.) in Table |^ 
using the alternative method of extraction described in the previous section. We 
remark that in this case (as one can check using Table I^J the inverse propagator is 
remarkably linear in k'^ up to the maximal (on-axis) momentum k'^ = 4. 

6 Conclusions 

In this paper we have shown that simulation data is in perfect agreement with the 
conventional renormalization group scenario in the lattice 1-component 01 model. 
This is contrary to the claim of ref. pQ. We have explained the deficiencies in 
their analysis. The quantitative agreement with analytically obtained results is 
completely satisfactory taking into account that some systematic errors are probably 
underestimated in ref. |3|. 

Conventional wisdom is not firmly established merely by a majority vote, but 
on solid theoretical and numerical studies. The same applies of course to deeper 
questions e.g. whether QCD is the correct theory of hadronic physics. The work 
reported here again reinforces the CW regarding triviality. Had the critique of the 
authors ^ been correct, they would have indicated serious non-standard implica- 
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tions concerning the Higgs' sector of the Standard Model. 
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Table 3: Measured values of x,^' from various Ising simulations; data from this investigation have 
no entry in the last column, kq = 0.080795. In all cases T = L except for the lattice where 
K = 0.076* denotes L = 20,T = 32. 
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16 


0.554(1) 


8.9 


0.563(1) 


0.886(3) 


36.0(2) 


0.076 


20 


0.392(1) 


7.8 


0.3950(9) 


0.896(4) 


30.37(28) 


0.076* 


20 


0.390(2) 


7.8 


0.3940(15) 


0.893(7) 


30.08(46) 


0.0754 


32 


0.266(3) 


8.5 


0.2666(25) 


0.932(18) 


27.1(1.0) 


0.0752 


36 


0.205(2) 


7.4 


0.2054(16) 


0.901(15) 


23.06(73) 


0.0752 


48 


0.205(2) 


9.8 


0.2055(18) 


0.905(17) 


23.21(83) 


0.0751 


48 


0.168(2) 


8.1 


0.1688(15) 


0.883(17) 


20.51(74) 



Table 4: Measured values of m, mn, Zr, gR for the Ising model. In all cases T = L except for the 
lattice where k = 0.076* denotes L = 20, T = 32. 
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0.3 


0.275376 


0.144499 


0.36895(2) 


0.5216(48) 


0.6 


1.177242 


0.130307 


0.13957(13) 


0.1873(15) 


1.0 


oo 


0.0751 


0.15657(4) 


0.1691(15) 


A 




9r 


X 




0.3 


0.962(12) 


20.2(4) 


12.40(5) 


0.975(20) 


0.6 


0.951(4) 


19.7(4) 


102.5(3.5) 


0.937(48) 


1.0 


0.896(8) 


20.6(7) 


206.4(1.2) 


0.886(17) 



Table 5: Results for simulations on the RG curve, gn ~ 20 
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